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Abstract 

We give an explícit stochastic Hamiltonian model of discontinuous 
unitary evolution for quantum spontaneous jumps like in a system of 
àtoms in quantum òptics, or in a system of quantum partides that 
interacts singular ly with "bubbles" which admit a continual counting 
observation. This model allows to watch a quantum trajectory in a 
photodetector or in a cloud chamber by spontaneous localisations of 
the momentums of the scattered photons or bubbles. Thus, the contin- 
uous reduction and spontaneous localization theory is obtained from 
a Hamiltonian singular interaction as a result of quantum filtering, 
i.e., a sequential time continuous conditioning of an input quantum 
process by the output measurement data. We show that in the case 
of indistinguishable partides or àtoms the a posteriori dynamics is 
mixing, giving rise to an irreversible Boltzmann-type reduction equa- 
tion. The latter coincides with the nonstochastic Schròdinger equation 
only in the mean field approximation, whereas the central limit yelds 
Gaussian mixing fluctuations described by a quantum state reduction 
equation of diffusive type. 
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Introduction 

This paper gives a microscopic theory of quantum jumps, inducing the spon- 
taneous localizations. It is based on a new type of quantum dynamics, de- 
scribed by the generalized (singular) stochastic Schròdinger equation, and on 
the quantum filtering process, amounting to a modification of quantum evo- 
lution along the quantum measurement trajectòries. We do not assume the 
existance of a unique universal mechanism for continuous reduction and spon- 
taneous localization during the measurement process. There are many such 
mechanisms, as many as the considered kinds of the singular interactions, 
amounting to completely unequivalent modifications of quantum evolution 
(e.g. quantum jumps and quantum diffusion) along lines, formally similar to 
different types of macro-objectification of pointer positions. 

The quantum measurement theory based on the ordinary von Neumann 
reduction postulate applies neither to instantaneous observations with con- 
tinuous spectra nor to continual (continuous in time) measurements. Al- 
though such phenomena can be described in the more general framework 
of Ludwig's Davies-Lewies operational approach [l]-[5] there is a particular 
interest in describing quantum measurements by concrete Hamiltonian mod- 
els from which the operational description can be derived by an averaging 
procedure. Perhaps the first model of such kind for instantaneous unsharp 
measurement of particle localization was given by von Neumann [6]. He 
considered the singular interaction Hamiltonian 

*,(«) = rf(t)^, íW = { ro ' (o.i) 

with Dirac 5-function, producing the translation operator 

s x = exp{~ j™ h x (t)dt} = e- xd ' AX (0.2) 

at time t — of the measurement. Here x is the position of the particle 
and A is the pointer position, defining, say, the momentum p of a quantum 
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meter. The particle scattering operator S = {s x } applied to the (generalized) 
position eigen-kets \x) as S'Ix) = \x)s x does not affect the position x of the 
particle but changes the meter momentum p to 

y = slps x = x + p, (0.3) 

which implies that the initial wave function ipo(x, A) of the system "particle 
plus meter" is transformed into 

ip(x, A) = s x ip (x, A) = if> (x, X-x). (0.4) 

If in the initial state the particle and the meter were not coherent, 
ipo(x,\) = r)(x)fo(\), and the wave function / of the meter was fixed, 
then one can obtain the unitary transformation S:^h s,ipo via a family 
{F(A)} of reduction transformations F(X) : i] ^ f·(X)r) for the prior particle 
vector-states r\. Specifically, Eq. (0.4) can be defined as 

i/>(x, A) = / (A - x)r)(x) = f x (\)r)(x) . (0.5) 

The linear nonunitary operators F(X) = s./o(A) act on \x) as the multipli- 
cation F(X)\x) = \x)f x (X) by f x (X) = / (A - x) = s x f (X) and would give a 
sharp localization of any particle wave function r](x) at the point x = X of 
the pointer position provided that the wave function /o(A) could be initially 
localized at A = 0. But there are no such sharply localized quantum states 
for the continuous pointer, and the best that one can do is to take an ap- 
proximately sharp wave packet /o(A). Say, / is the Gaussian packet rescaled 
to the Standard form 

/o(A)=exp{-|A 2 }. (0.6) 
This results in the unsharp localization [F(X)rj](x) = ip(x, A), 

^AHexpj-^Or-A) 2 }^) (0.7) 

of any wave function rj(x) about the observed value A of the pointer. The 
localized wave function if)(x, X) defines for each measurement result A the 
posterior particle vector-state 

ifc = cíV(A), |ca| 2 = ||^(A)|| 2 (0.8) 
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up to the normalization 



||^(A)|| 2 = J |^(x,A)| 2 dx = p(A) 



(0.9) 



to the probability density 



p(A) = I \f (\-x)\ 2 \ V (x)\ 2 dx 



(0.10) 



For commuting operators x and y, this is equivalent to the classical mea- 
surement model y = x + p for unsharp measurement of an unknown signal 
x via the sharp measurement of the signal plus Gaussian noise p with given 
probability density |/o(A)| 2 of the random vàlues p — X. 

These arguments illustrate how to intèrpret the reduction model involv- 
ing the continuous spectrum of a quantum measurement as a Hamiltonian 
interaction model with nondemolition observation for a quantum object us- 
ing the measurement of the pointer coordinate of the quantum meter. Since 
von Neumann introduced this approach, they were used by numerous other 
authors [7, 8] for the derivation of a generalized reduction r\ \— > F(X)r) that 
would replace the von Neumann postulate rj \— > E(\)r) given by orthoprojec- 
tions {E(X)} in the case of discrete vàlues A. 

As extended to nondemolition observations continual in time [9]-[15], this 
idea constitute an essence of the quantum filtering method for the deriva- 
tion of nonlinear stochastic wave equations describing the quantum dynamics 
under the observation. Since a particular type of such equations has been 
taken as a postulate in the phenomenological theory of permanent reduc- 
tion, quantum jumps and spontaneous localization [16]-[20], the question 
arises whether it is possible to obtain this equation from an appropriate 
Schròdinger equation. Here we shall show how this can be done by second 
quantisation of the singular interaction Hamiltonian considered by von Neu- 
mann obtaining a stochastic model of continual nondemolition observation 
for the position of a quantum particle by counting some other quanta. 

First we show in the Sec.l that even the projection postulate can be 
derived in the framework of this approach with the suggested Hamiltonian 
interaction and a proper nondemolition observation with a discrete spectrum 
and sharp initial state of the meter. Then we define a single-kick stochastic 
wave equation for the reduced state-vector, corresponding to the unsharp 



/o(A). 
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In the Sec. 2 we develop the differential treatment of discontinuous uni- 
tary evolution in terms of generalized Schròdinger equation corresponding 
to the scatterings at given or randomly distributed time instants. Then we 
show how it can be reduced to the many-kick stochastic wave equation, de- 
scribing spontaneous localization of a quantum particle under the continual 
observation of its trajectory in a bubble chamber. 

The quantum system of many similar interacting partides in a bubble 
chamber is treated in the Sec. 3. We prové that the reduced dynamics of 
the partides is mixing under the continual observation of positions of the 
scattered bubbles. It can not be described by any reduction equation for the 
posterior wave function but can be described by an irreversible stochastic 
equation for the density matrix of the partides. 

Finally we obtain in the Sec. 4 the macroscopic and diffusion limits of 
the generalized Schròdinger and reduction equations under the limit of weak 
coupling constant and higher frequency of the interaction with the appara- 
tus. In the macroscopic limit these two equations coincide with the ordinary 
Schròdinger equation, corresponding to the mean field approximation, while 
in the diffusion approximation, which describes the fluctuations, they lead to 
the essentially different kind of stochastic wave equations [23, 24]. 

These results are presented by the equations (1.10), (2.10), (3.10) and 
(4.10). 

1 Hamiltonian Model for a Generalized Re- 
duction 

Let TL be a Hilbert space called the state space of a particle, and let R be a 
selfadjoint operator in TL with either integer or continuous spectrum Z or R. 
Let k > be a scaling parameter. One can regard the scaled operator kR as 
the position x of the particle in R or in the lattice kZ if it is quantised so that 
kR is given in the x-representation, the operator acts as the multiplication 
kR\x) = \x)X(x) by \(x) = k[x/k\, where |_^J £ Z denotes the integer part 
of x. 

A quantum meter with continuous (A = R) or discrete (A = eZ) pointer 
scale is described by the Hilbert space L 2 (A) of complex- valued functions 
/ : A — > C square integrable in the sense that ||/|| 2 = J\f(y)\ 2 dy < oo. 
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The last integral coicides with £ yeA \f(y)\ 2 e if f(y) = T,iez8f(y)f(el)£ is 
the isometric interpolation of the discrete wave-function f(el), l G Z, given 

by 

Sf(y) = f 1/£ e 27TÍk(£l - y) dk, yeR. 
Jo 

Consider a moving particle with Hamiltonian H in 7i. Its singular evo- 
lution corresponding to the position measurement at time t = is described 
in the product space 7íi = 7í® L 2 (A) by time-dependent Hamiltonian 

Hi(t) = H ®1 — kR® S(t)Q , (1.1) 

in the interaction picture corresponding to the free translation evolution of 
the poiter coordinate pGA. Here 1 is the identity operator in L 2 (A), defining 
the particle Hamiltonian H = H(B) 1 in ?í 1 , and Q is treated as a coordinate 
operator of the meter, Q = ihd/dy, having the bound spectrum 2nh[0, 1/e) 
for A = eZ. This generates the shift operator (0.2) in L 2 (A) with x G A and 
the scattering operator 

* = «p(í«r®«,) = {?; HI d.2) 

in Tïi, where Q t = ^tQ, 1* = 1 if í > 0, and l t = otherwise. 

The singular time dependence of the Hamiltonian (1.1) makes it impos- 
sible to define the Schròdinger equation ihdip/dt = Hi{t)ip{t) at t — 
in the usual sense. But one can define a discontinuous unitary evolution 
Ui(t) : Tii — > Tíi for some t < as the single-jump unitary process 

C7i(í) = exp(-^ J* H^ds) = e iH ^ t0 ~ t)/n S t 

provided that [R, H] =0. If i? and H do not commute, the usual method 
which a physicists would follow at facing this difliculty is to consider the free 
evolution ip-(t) for t < and ip+(t) for t > 0, supplemented by the boundary 
condition tp+(0) = Sip-(0). Although this method is hardly applicable in 
the case of spontaneous interaction at a random instant íi > tç,, let us show 
that the discontinuous unitary evolution corresponding to the fixed íi = 
can be defined as the fundamental solution ^(í) = Ui(t)ipo of a generalized 
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Schròdinger equation. Such regularized equation can be written in terms of 
the forward differentials 

d^(í) = V(í + dí) - V(í), dl t = U+dt - U , 
as the Ito integral of 

dij;(t) + ^H iP(t)dt = (S-I)iP(t)dl t , ip(t )=^ . (1.3) 

For every initial ipo e T~í the equation (1.3) has a unique solution at ío < 0; 
this solution is given by the unitary operator Ui(t) = U (t —t ) S t (—t ) , where 
U (t) = exp{-iH t/h}, S t (r) = U^(r)S t U (r). 

To prové this let us rewrite the generalized Schròdinger equation in the 
integral form 

^(í) = e'^ot/h ^H t /h^ + J* e iH r/h( S _ J)^( r ) dl ^ ( L4 ) 

(the equivalence of (1.4) to (1.3) can be shown by straightforward differen- 
tiation). We can write ip(0) = Uq{— t )ip for the solution (1.4) when t — 
since J t ° (p(r)dl r = if t < 0; hence ip(t) can be rewritten as 

iP(t) = U o (t)(iP(0) + (S- 7)1^(0)) = U (t)SM0) . 

It follows from / t * </?(r)dl r = 1^(0) for any function </? and any t < 0, and 
(S — I)l t = S t — I by the definition of the scattering operator (1.2). Thus, 
Ui(t) is the unitary operator 

U (t)S t U (-t ) = U (t - t )S t (-t ). 

This gives the solution to equation (1.3) for t = as well, since if;(t) = 
U o (t)S t ip(0) is equal to ipo for t — 0. 

The rescaled pointer operator P = k~ 1 I ® p, switched on at the instant 
t — of the scattering, is described in the space H <g> £ 2 (A) by the operator 
Pt = UP, having in the Heisenberg picture Y t = U\{t)P t Ui{t) the constant 
form Y t = Y(-t ) for all t > 0. Here Y(r) = £fi(r)PS(r) taken at r = -í„, 
t < is the shifted operator 

Y(r) = R(r)®l + -I®p, (1.5) 
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where R(r) = U'(r)RU(r) and S(r) = exp(KR(r)<giiQ/h). Obviously the ob- 
servables {Y t } are self-nondemolition in the sense of their joint measurability. 
It follows from the trivial commutativity condition 

[Y s ,Y t ]=0, Vs,í. (1.6) 

And they are nondemolition with respect to an arbitrary particle operator 
X t {r) = Ui(ty(X <S> l)Ui(t) in the Heisenberg picture at t — ~ r m the sense 
of their predictability [9, 13] 

[X s ,Y t ]=0, Vs>t, (1.7) 

indeed, [X 8 , Y s \ = and Y s = Y t for s, t > since 

UlWPMs) = &(-t )P t S(-t ), Vs > t. 

Let us fix a state vector fo G L 2 (A), \\fo\\ = 1, given by a wave function 
fo(y) on A localized at y — 0. Let |y) denote the (generalized) eigen func- 
tion in the spectral representation p = J y\y)(y\dX of the pointer coordinate 
(pf)(y) — yf(y)i where dÀ = dy if y G R and dÀ = e if y G eZ. This yields 
the localizing transformations F t (y) = (y\S t f in the form 

«<») = /.<«" = ^° (1.8) 

where F(y) = /o(y^ — t_R). The reduction transformations r] \— > ip(t,y), 
defined on the particle space by the formula 

ïP(t,y) = U(t-t )F t (-t ,y)ri, y G A, 

with C/(í) = exp(—iHt/H) and F t (r) = W{r)F t U{r) reproduce the unitary 
evolution Ui(t) in the y-representation on r\ ® f G fíi similarly to Eq. (0.4) 
and (0.5), namely, 

U(t - t )F t (-t Q , y)v = U(t)(y\S t f U(-t )r) = (y|£/i(í)(?7 ® /o) ■ 

The operator F(r,y) = fo(yI — nR(r)) at r = — 1 with a given initial 
vector-state rj of the particle before the scattering (í < 0) define the proba- 
bility measure 

»(A)= f ||F(r, í /)7 ? || 2 dA = (r ? ,n í (A)7 ? ), AÇ A, 

JA 
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for the statistics of the nondemolition measurement of kR via the observation 
of the pointer position y G A after the scattering. It is given by a positive 
operator-valued measure II t (A) = J A \F t (y, —t )\ 2 d\ at t > 0. (Before the 
scattering it defines the initial probability measure 

Aí (A) = (j7,n t (A)77)= / \fo(y)\ 2 d\, t<0 

JA 

which is independent of 77.) 

In the case of discrete scale the eigen functions \y) are normalisable: 
(y\y) = an d one can take a sharply localized /o = e d//2 |0), given by the 
eigen function fo(y) = e(y)/e d ^ 2 , where e(y) — 1 if y — 0, e(y) = if y 7^ 0. 
By renormalising the operator F(y) as E'(y) = e d ^ 2 F(y) with step e = k of 
the position quantisation, one obtain the orthogonal projections 

e(nR — yl) = \ \x)(x\dx = E(y) , y G k'L 

J x:\_x / n\=y / k 

of kR = J \(x)\x)(x\dx = J2 y e K zyE(y) corresponding to eigen vàlues X(x) = 
k[x/k\. Thus, II t (A), A Ç A is the spectral measure J2 y eA E(to, y) of the 
quantised position KR(-to) of the particle for a t > given by the eigen or- 
thoprojections E(r,y) = W(r)E(y)U(r) of the operator R(r), corresponding 
to the rescaled pointer integer vàlues y/K. Thus, the projection reduction 
postulate has been deduced from the Hamiltonian interaction (1.1) and the 
nondemolition measurement for the sharply localized initial state /q. But 
there is no continuous limit as k — > of such sharp reduction with non- 
trivial e ^ 0, since ||e|| 2 = j \e(y)\ 2 d\ — k — > and the sharp function 
e(y) = 5Q(y)n on A = kZ disappears as an element of the state space £ 2 (R) 
of the continuous meter. 

In the continuous case, one can take only an approximately sharp f G 
L 2 (R), say fo(y) = 5g(ï/)e, y G R, and renormalise the operator (1.8) as 
Gt(y) = fo(yI - K Rt)l 'fo(y) if fo(y) °, as in the Gaussian case (0.6). They 
define G t (y) as the identity operator for t < or r > 0, whereas for t > 
and r < the operator 

G(y) = (y\Sf = (y\G, G = f 1 Sf , (1.9) 

say, of the Gaussian form 

G(y) = exp{7cnR(yI -^kR)} 
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given by the generalized eigen functions \y) = \y)/fo(y) for the spectral rep- 
resentation p = J y\y)(y\dfi with respect to the initial probability measure 
d/zo = \f (y)\ 2 d\ with density \fo(y)\ 2 = exp{— ny 2 } in the Gaussian case. 
The corresponding propagators 

T(t, y) = U(t - to)G t (-t , y), yeA,t <0, 

where G t (r) = W(r)G t U(r), define the operator-valued measure II t (A) as 

n t (A)= / T\t,y)T(t,y)à^= í \G t (-t ,y)\ 2 dfi 

JA JA 

so that the output probability measure \i is absolutely continuous with re- 
spect to fi for any 77 G 7í: //(A) = J A \\G(— 1 , y)r]\\ 2 dfi . Hence, the reduced 
state vector x(í, y) = T(t : y)rj is normalized to 1 as a stochastic vector pro- 
cess x(t) '■ U l— ¥ x{t->y) G ^ in the mean square sense with respect to the 
input probability measure /io, 

\\x(t)\\l = J (x(í,y),x(í,y))d/x = (77,77) = 1 . 

This model of nondemolition observation with continuous data y G R also 
applies to unsharp measurement of operator R with discrete spectrum. In 
contrast to sharp measurement, unsharp measurement is not sensitive to the 
continuous spectrum limit as k — > of R = [x/k\, corresponding to the 
replacement of kR by x G R. 

For any initial t ü <Q the stochastic vector process x(t) = T(t)rj satisfies 
the single-kick equation 

d X (t) + ^H X (t)dt = dl t [G-I]x(t), X (to)=V (l·lO) 

generated by the random differential dlf[G — I\{y) = (G(y) — í)dl t on y G A 
with respect to the initial probability measure fj, . This simplest reduction 
equation is written in terms of the forward differentials d%(í, y) = x(t + 
dí, y) — x{t-iy)-> that is, is understood in the sense of Ito. 

Indeed, by representing Gt as Gt{y) = (G(y) — Í)lt + I = I + lt[G — I](y) , 
one can describe X (t) —U(t — t )G t (—to)rj by the integral equation 

X (t) = U(t - t )v + U(t)l t [G - I]U(-t ) v = 

= e- iHt ' h (e™*'\ + £e™ r ' R dl r [G-I] X (rj) . 
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But this equation is equivalent to the differential equation (1.10), a fact that 
can be proved by straightforward differentiation taking into account the Ito 
multiplication table 

(dt) 2 = 0, dídl t = = dl t dí , (dl t ) 2 = dl t . 

Similarly, one can the simplest nonlinear stochastic equation for the nor- 
malized reduced state vector X y(t) = x(ti y)/\\x(ti 

d Xv (t) + ^H Xy {t)dt = (G y (t) - I)xy(t)dl t , x v (0) = V 

where G y (t) = G(y)/\\G(y) Xy (t)\\, rj e H. This equation is an equivalent 
differential form the nonlinear integral stochastic equation 

Xy (t) = e - m /*(^\ + jíV··/ ft (G w (r)-/)x w (r)dl r ) = 

= U(t - t )rj + U(t){G y - I)hXy(0) = U(t)G t (y) X y(0)/\\G t (y)xy(0)\\ . 

This yields X y(t) = T(t,y)r}/\\T(t,y)r}\\ for aí o <0 because of \\G t {y)r]\\ = 
\\nt,y)ril 

Note that the random state vector Xy(t) ^ s obtained by conditioning 
with respect to the output (rather than input) probability measure d/i = 
llx(í,í/)l| 2 d/io. 

2 Spontaneous Localization of a Single Par- 
ticle 

Let us consider a spontaneous process of scattering interactions (1.1) of a 
quantum particle (or an atom) at random time instants t n > 0, t\ < t 2 < . . ., 
with a renewable meter in an apparatus of the cloud chamber (photodetector) 
type with bubbles (photons) serving as the meter. We consider the increasing 
sequences (íi,Í2, • • •) as countable subsets r C R + such that r t = r fi [0, t) 
is finite for any t > in accordance with the finiteness of the number of 
scattered bubbles on the finite observation interval [0, í). The set of all such 
infinite r will be denoted by T^, and T is the inductive limit UT t as t — > oo of 
r t = {r t \r G Too}, which is the disjoint union r t = S^L o r t (n) of n-simplice 
r t (n) = {íi < ... <t n } C [0,í) n . 
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The measurement apparatus is assumed to be a quantum system of in- 
finitely many bubbles each of which is identical to the single meter described 
in the previous section. The pointer coordinate is attached to the momentum 
p n of a bubble labeled by the scattering number n G N, such that it shows 
the momentum À n G A of the scattered bubble at each time t n G r. 

The corresponding Hamiltonian of the moving particle is given by the 
series 

oo 

H(t,T) = H -KR®J2 S ( t - t n)Q(n) (2-1) 

n=l 

having at most two nonzero terms if t G r. Here H = 1®°° is the Hamil- 
tonian describing the time evolution on the intervals between the scatterings 
í G r and Q(n) is the coordinate of the nth scattered bubble, given as the 
operator Q(n) = ihd/d\ n . 

The generalized Schròdinger equation corresponding to the Hamiltonian 
(1.1) can be written for fixed r G by analogy with the single-kick case 

dV>(*) + \ H il>(t)dt = (S(n t ) - I)ip(t)dnt, V(0, r) = ■ (2.2) 

Here S(n) = exp{^/ti?® Q(n)} and n t {t) = |r t | is the numerical process that 
gives the cardinality \r t \ — Y^rer h-r of the localized subset r t = {t n < t}, so 
that dn t (r) is equal to 1 for t G r, and zero otherwise. 

The solution to this equation is uniquely determined for every r G by 
the initial state ipo of the system. Namely, 

Íj(t,r) = U(t,r)^ Q , 

where U(t,r) = U (t)V t \T), V/(r) is the chronological product YiTer $t(r) = 
S(t nt ) . . ■-S'(íi), and 

Vi(r)=5Í(í 1 )5Í(í 2 )...= (n5 t (r)) t . (2.3) 

Ver ' 

Here S t (t n ) = Ul{t n )St{n)UQ{t n ) for t n < t, where S(n) = exp{ — ^ kR ® 
Q(n)}, and ^(ín) = / if t n > t so that the infinite product (2.3) contains 
only a finite number n t = J2 r er ^t-r of factors different from the identity 
operator /. 
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Recali that the differential equation (2.2) is equivalent to the integral 
equation given by the recurrence relation 



^(í, r) = e- iH ^ n U + g ^ h (S(nr) - /)^(r, r)) (2.4) 



for every r G T^. Hence, ^(í,r) = U (t)Vf(T)if> , where Z7 (í) = e _ifl °*/ R 
and V t (r) is a solution to the operator equation 



W = / + £ K(í)(4(r, r) - /) , ^o(r) = / , 



where S^í, t) = Uo(t)^ S(nt(r)yUo(t). But this equation has a unique solu- 
tion (2.2), which can be written as the binomial sum 

[L t (h) + I] [L t (t 2 ) +/]...= £ L(si, r) . . . L(s n , r) 

o-Çtí 

in terms of er = {si, . . . , s n }, si < . . . < s n , n < n t , L t (r) = Sj(r) — I (= 
if r > t) and L(r, r) = S^r, r) — /. Indeed, this sum contains / as the null 
product corresponding to er = 0, and the sum of the other terms is equal to 

r<t 

V t {r)-I = J2T, L (si,r)...L(s m ,r)L(r,r) = 

rer crÇTr 

= g V r (r)L(r, r) = g VÇ(r)(fift(r, r) - J) , 

where m < n t — 1. 

Note that the generalized differential equation (2.2) depending on r e 
via n t = n t {j) is not necessarily stochastic as long as we have not fixed 
a probability distribution for the instants r = (ti,t 2 ,···) of the singular 
interactions. In order to obtain a continuous (at least in the mean) dynamics 
for such a quantum jump process it is necessary to assume the interactions are 
spontaneous with a continuous probability distribution of random instants 
t. One can assume that the number process n t (r) is stochastic, given by the 
Poisson law 7r (dr) on presented as the projective limit as t — > oo of the 
probability measures 



7T 



(dr t ) = e'^u^dn, u>0. (2.5) 
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Here r< = r is a finite time-ordered sequence r(n) = (íi, . . . ,t n ) G T t with 
n = n t , dr t = YikL·i is the measure on T t given by the sum of product 
measures dti, . . . ,dt n = dr(n) on the simplices T t (n), dr(0) = 1 on r t (0) = 
{0} such that 

í v ^dr:= V/ í ... / díi . . . dt n = e ut . 

Jv t n=0 J Jo<t 1 <...<t n <t 

Note that any other numerical process can be described by a positive density 
function /(r) with respect to the Poissonian measure, that is, has the form 
f(r)n(dr). 

Let us fix an initial state ipo = /q° of the bubbles as the infinite product 
f£° = <S>kLifk of the identical state vectors f\ = fo of the bubbles given by 
a normalized element / G L 2 (A). This defines the solutions ip(t,r) of the 
stochastic equation (2.2) with the initial data ip = r\ <g> </? given by state 
vectors i] £ ?í in the particle space as the state vectors in the product 
space Tíoo = 7í <g> £, where £ is the Hilbert space generated by the infinite- 
product functions <f(v) = O^Li /fc(^fc) with the equal elements fk = fo for 
almost all k. We suppose, as in Sec.l, that ipo(v) ^ for almost all v = 
(Ai, À 2 , . . .) such that one can identify the space £ with the space Lo(A°°) of 
all square-integrable functions / = y?/Vo with respect to the product measure 
/x^°(dA) = /i (dÀi) •/!o(dÀ 2 ) • • • on the space A°° = A x A x . . . of the sequences 
A= (À!,A 2 ,...): ||/||§ = /|/(A)| 2 ^(dA) = Hv^ll 2 . The generalized product- 
vectors |A) = |Ai) <g> |A 2 ) <E> • • • of this space are defined by tensor-product of 
the 5-functions (A'|A), respectively to /i such that / | A) (A| / ug°(dA) = l® 00 . 

Consider the sequence (pi,p 2 ,...) of momentums p n , of the scattered 
bubbles at the time instants {íi,í 2 , • • •}• The commuting operators p n , n G 
N, described in £ by the multiplications p n \\) = |A)À n , are assumed to be 
measured at the random time instants t n , n G N. The point trajectòries 
of such measurements are given by the sequences = (yi,y 2 ,...) of pairs 
Un = (t n , A n ) with ti < t 2 , . . . and A n G A, identified with countable subsets 
{yi, y 2 , . . .} C R+ x A. As elements v = (r, A) of the Cartesian product Yqo = 
Tqo x A°°, they have the probability distribution P (dv) = 7r (dr)/ig°(dA), 
where A°° is the space of all sequences A G A°°, À n G A, equipped with the 
probability product-measure /xg°(dA). 

The measurement data of the observable process up to a given time in- 
stant t > is described by a finite sequence v t — (yi, . . . , y n ) with n = n t (v) 
given by the numerical process n t (r) for the component r off. 
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Let us introduce the counting distribution n*(A) = l^n (R+ x A)| as the 
number n t (A,v) of scatterings in the time-space region [0, t) x A and define 
the counting integral J °° J A L(r, À)dn r (dÀ) over y G R + x A as the series 

n[L]{v) = Y,L(y), V.eL. (2.6) 

y£v 

Having fixed an integer-valued distribution n t (A) G {0, 1, . . .} as a function 
of t > and of measurable sets A Ç A, one can obtain the corresponding 
trajectory v as a sequence of the counts of the jumps of n t (A) in the time- 
space R + x A. 

Given an initial state vector in Hoo of the form ijj = r\ <g> <p with fixed 
V^o = /o°, one can define a nonunitary stochastic evolution r] \— > T(t,v)r) by 
setting 

T(í,t;) = (A|f/(í,r)v9o, u=(r,A), 

which reproduces the unitary evolution C/(í, r) = Uo(t)Vt(r) defined by (2.2). 
This can also be written as T(t,v) = U{t)F}(v), since U (t) = U(t) <g> 1®°° 
commutes with the (generalized) eigen-bras (A| = (Ai, A2, . . . | of the bubble 
coordinates (pi,P2, • • •) : (^l^o(^) = U(t){\\. The reduction transformations 
F t (v), v G Tqo, are given by the chronological products 

F t (v)=G\( yi )G\{y 2 )... = Y[G\(y) (2.7) 

yev 

of G t (í n ,A n ) = W(t n )G(\ n )U(t n ) for t n < t, where G(X) = (X\Sf , owing 
to the product form (2.3) of the unitary transformations V t (r), r G T, and 
(A| = ®r=i(A fe |, <PoW = nr=i /o(A fc ), (A|^> = l for A G A°°. The stochastic 
operator (2.7) defined by the single-point reductions 

G t (r,A) = (A|5 t (r)/ = | ^ A) ' (2.8) 

is normalized with respect to the initial probability P (dv), n Too [/](í) = /, 
where 

Il A [X](t)= í T*(t,v)XT(t,v)P (dv)= í F t (v)XF}(v)P (dv), (2.9) 

JA JA 

is a continual operational-valued measure [3]- [5] defined on measurable sets 
^4 Ç Tqo of the point trajectòries Vt = {(r, A) G i>|r < £} given by the 
operations 3>t(i>) : X 1— > F t (v)XF/(v) for particle operators X : — > 7i. 
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The positive operator-valued measure II t (A) = ITi[.f](í) gives the statis- 

tics 

P t (A)= í \\(X\U(t,r)i, \\ 2 ^(àX)n (dr t ) 
JA 

for the continual observation with respect to an arbitrary initial wave- 
function ip — 7] <S> /o°, r] G TC in the form 

p t (A) = (77, n t (A)77) . 

The output probability measure P(A), A Ç Yoo, is defined by the marginales 
P t (A), A Ç T t , as t -> 00. 

The reduced wave function xí^^) — T(t,v)r] is normalized 

IIX(0II 2 = / ||x(í^)|| 2 Po(d^) = l 

as a stochastic vector process x(í) : — > 7i with respect to the initial 
probability distribution P on T^. 

It satisfies the stochastic wave equation 

d X (t) + ^ #x(t)dí = dnt[G - I] X (t) , x(0) = V , (2.10) 

expressed in terms of the random differential dn t [G — I](v) = (G(X nt (v)) — 
l)dn t (v), n t (v) = n t {r) for the point distribution n t [L] = J A L(X)n t (dX) = 
n[L t ] over À G A defined in (2.6) for L t (r, X) = l í (r)L(À). 

To prové Eq.(2.10), discovered for the first time in [13], we rewrite it in 
the following integral form 

X(t) = e-™'*^ + fjj Hr l\G{X) - J) X (r)dn r (dA)), 
given for each v G by the finite sum 

n[l t rf(r)(G-l) X ](v)= £ ^(r)(G(X) - I) x (r) . 

(r,X)ev 

We express the solution to this equation in the form x{t> v ) — U(t)F?(v)i] 
via the solution (2.7) to the recursion equation 

F t (v)=I+ ]T F r (v)(G(r,X)-I), F (v)=I, 

(r,\)ev 
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with G(t, X) = U'(t)G(X)U(t), as was done for the unitary case. 
Let us also write on F the nonlinear equation 

d Xv (t) + l - H Xv (t)dt = (G v (t) - I) Xv (t)dn t (v) , Xv {0) = V, 

with G v (t) = G(\ nt ( v \)/\\G(\ nt r v ))Xv(t)\\- Its solutions define the normalized 
reduction Xv (t) = X (t,v) /\\x(t,v)\\ for continual counting measurements as 
a stochastic vector process Xv {t) G 7í with respect to the output probability 
measure P of the point process t \— > v t . This can easily be obtained, as in 
[14] by applying the Ito multiplication table 

(dí) 2 = , àtàn t = = dn t dt , (dn t ) 2 = dn t . 

3 Mixing Reduction for Many Identical Par- 
tides 

We now consider M identical partides (àtoms) interacting independently 
with the bubbles (photons) in accordance with the scattering term in the 
Hamiltonian (2.1). The spontaneous process of scatterings is described by the 
time-ordered sequences of pairs (k n , t n ), ti < Í2 < • ■ •, where k n G {1, . . . , M} 
is the number of the particle labeled by the scattering number n G N at time 
instant t n > 0. We have excluded possibility of two or more scatterings of the 
bubbles at the same instant of time, as was done for a single particle in Sec. 2. 
The sequence (ki, íi), (k 2 , t 2 ), ■ ■ ■ of the scatterings can be represented by the 
occupational subsets = {t n G r\k n = k} of the time set r = {ti,t 2 , . . .}, 
which are disjoint, Tk H T\ = if k ^ l since the scatterings for different 
partides are independent. We shall consider the M-tuples r. = (ti, . . . , tm) 
of these countable subsets r k C R+ as elements r. G of the Cartesian M- 
product of Toc, given by the partition r = Ut^ : = Ur^, fi T\ — if k ^ l of 
a t G Toc. 

The Hamiltonian of the interacting partides for fixed r. = {(fci,íi), 
(k 2 ,t 2 ), . . .} reads 

oo 

//(í, r.) = F M - k 53 i?(fcn)) ® 5(t - í n )Q(n) . (3.1) 

n=l 
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Here Hq = H ® 1 is the Hamiltonian of the partides describing the time 
evolution on the intervals between the scatterings with the bubbles: 

M MM 
k=l k=ll>k 

where H(k) = J®^" 1 ) <g> H <g> I®( M - k ) i s the Hamiltonian of the fcth particle 
and W{k, ï) is the interaction potential in H® M of the fcth and Zth particle, 
1 < k < l < M. 

Let 7í M denote the M-particle Hilbert space, which is an invariant sub- 
space 7i M Ç 7i® M of symmetric (bosons) or antisymmetric (fermions) M- 
tensors r/ M G Tl® M generated by the product-vectors ®^ =1 rjk G Tl® M with 
r/k G H. The correspondent Ito-Schròdinger equation for the stochastic state 
vector ip M (t) : r. i— > ip{t,r,) with vàlues ip(t,T,) G in the product space 
= H M <g> £ of the M-particle space 7Í M by £ = lim^*, L 2 (K n ) reads 

d^ M (t) + i ií M ^ M (í)dt = (S-ífct, n t ) - I)^ M (t)dn t . (3.2) 

Here S(k, n) = exp{i rei2(Jfe) <g> Q(n)}, R(k) = J^" 1 ) <g> i? ® /®(Af-fc) > 

= E , wh ere l T (í) = | | í £ r 

is the random number fc t : — > {1, . . . , M}, labeling a particle by at any 
instant t G r k of its collision with a bubble labeled by ra^r) = Y^k=i n k,t — 
\rD [0,í)|, where n M = |r fc n [0,í)|, r = Ur fe . 

The solutions ■?/>(£, r.) = U(t, r.)^^, G 7Y^, of equation (3.2) can be 
written as C/ (í, r.) = [/^(í)^ ( T ·) m terms of the finite chronological product 

V t (r.) = Sl(k 1 ,t 1 )Sl(k 2 ,t 2 )..., r.GC (3.3) 

where S t (k,t n ) = Iiït n >t, S t (k,t n ) = U™\t n )S(k : n)U™ (t n ) if t n < t, and 
U^(t) = exp{ — -i- H^(t)}. The proof is exactly the same as for the case of a 
single particle M — 1. 

Let cu = (u>i,u>2, • • •) denote a chronologically ordered sequence of triples 
w n = (k n , t n , X n ) and íl the space of all such sequences with {íi, t%, . . .} G r^. 
Every sequence uj G Í2 can be represented as a pair cu = (r., A), where r. = 



18 



(ri, . . . , tm) is a partition of the corresponding sequence r = {íi, t 2 , • • •} and 
A = (Ai, A2, . . .)) so that Q can be identified with the product x Too. The 
space Q is equipped with the probability measure P (du;) = 7ro(dr.)/ig°(dA), 
where 7r (dri, . . . , (Itm) = Ilfcli ^o(^- T k) is the product of the identical Pois- 
son measures (2.5), in accordance with the independence of the spontaneous 
interactions of each particle with the bubbles. 

Given an initial state vector ip M = i] M <g> ip , where ip = f£°, one can 
easily prové that the nonunitary stochastic evolution 



T(t,u) = U M (t)F?(cü) , F t (u) := G\{ Wl )G\{w 2 ) .... 

Here U M (t) = exp{-jrH M t}, G t (k,t n ,X) = I for t n > t, and G t (k,t n ,X) = 
U M ^ (t n )G(k, X)U M (t n ), t n < t, is defined by the reduced scattering operator 

G(k, X) = J®**" 1 ) ® G(X) ® I^ M ~ k) , G(X) = (X\Sf (3.4) 



for /o G L 2 (A), || /o|| = 1) applied to the A;th particle only in 7í M . The 
stochastic operator T(t) defines the solutions x M (t,iu) = T(t,u)r] M to the 
Ito differential equation 



for the stochastic vector states x M if) '■ ^ ~~ ^ M °f the M-particle system, 
corresponding to an initial f] M G 7i M . The right-hand si de of this equation 
is written as the forward increment of the point integral 



given by the stochastic distribution n[L\{uj) = ^2 w€u) L(w) for L t (k,r,X) = 
l t (r)L(k, r, X). The vector x M (t^) e ^ M as well as x(t,t>) in (2.10), is no 
longer normalized (||x M (í,u;)|| 7^ 1) for an initial state-vector i] M G 7í M , 
\\t] m \\ = 1, but it is normalized with respect to the probability measure Po 
on Q in the mean square sense. But, in contrast to x(ti v )i X M ifi ÜJ ) líi n °t 



T{t,u) = (\\U(t,r.) Vo , u = (t.,\) 



is also a finite chronological product 




x M (o)= v 



M 
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yet the reduced description of the M-particle system under the observation 
of the scattering process v t = {(r, À) G v\r < t}, given by the registration of 
the pointer positions \ n G A at random time instants t n . 

The reduced dynamics corresponding to the observation is described by 
a stochastic operational process X i— > 9[X](í), 

Q[X](t,v) = <S> t [U M \t)XU M (t)](v) (3-5) 

for the M-particle operators X : 7í M — > 7i M given by the conditional expec- 
tat ion 

$ í [X](r,A) = -l ïï £ F t (a.,X)XF t \a.,X) (3.6) 

Ucr fc =r 

where the sum is taken over all partitions <r, = (oi, . . . , o'm) of a finite subset 
T| = r fi [0, í). This averaging is due to the impossibility to detect the 
individuality of the identical partides producing the indistinguishable effects 
on the bubbles by measuring the scatterings of the bubbles. 

To prové Eq. (3.6), we need to compare the correlations of F t (u)XFf(u) 
and of an arbitrary functional g(v t ) of the observable point process v t with 
the correlations of (3.6) and of g(v t ). But by applying the well known formula 
[21] one can easily find 

M 

lx(a 1 ,...,a M )Y[da k = x(oi, . . . , o M )do 

JV t k=l Vt Vk--Lla k =a 

for the múltiple point integration that these correlations with respect to the 
probability measure Po on Q given by the Poisson law (2.5) simply coincide: 

,7r t (do.) / g(utí 1 a k ,X)F(a.,X)XF t \a.,X) l· i^(à\) 

„ M 

= / <í?(uf =1 o fe ), X{a u o M ))o II e-^^da k 

I T'M 

JL t k=i 
= í (g(a), Y X(o 1 ,... ) o M ))oe^VHdo 

4 crfc:Ucr fe =(T 

= f r ^(àff)f A _g(a,\)-^ E F(u)XF\u;)^(dX). 



Here (., .)o is the abbreviation for the inner product in £ of the test func- 
tion v i— > g(r,X) with fixed r G T t and the operator function X t (r,,X) = 
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F t (r„ A)IF/(t„ A) with fixed r. G rf. The probability measure 

M 

<(dr t )= £ U^(da k )=e- M ^\Mu\^dr t (3.7) 

Ua k =T t k=l 

on r t has the intensity Mv. It is induced by the measure 7To(dr.) on with 
respect to the particle identification map r. G T^f i— > r = Ufc=iTfc, defining 
the observable data v = (r, A) by the stochastic map uj = (r., A) i— > v G Too 
on G íl 

By the coincidence of the correlations proved above, the stochastic oper- 
ator (3.5) is indeed the conditional expectation of the stochastic operators 
X t {ui) with respect to the observable process v t . 

In contrast to the pure operations X \— > X t (u>) the reduction operation 
X i— > $ t [X](v), preserves the symmetry of the M-particle operators X with 
respect to particle permutation. It is the least mixing operation which pre- 
serves the indistinguishability of the partides respectively to the observations 
of the bubble scatterings, corresponding to the complete nondemolition mea- 
surement of the partides. 

Indeed, the reduction operation (3.6) can be simply written as the finite 
iteration 

1 M 

jj- G t(h,yi)---Gl(k n ,y n )XG t (k n ,y n )...G t (k 1 ,y 1 ) 

= %[. . . %[X]( Vl ) . . ]{y n ) = <S> t [X]( yi , ...,y n ) 
with n — \r t \ single mixing reductions 

i M 

nX](y) = -J2G\k,y)XG(k,y). (3.8) 
m k=i 

Given as the arithmetric mean value of the permutations for the pure opera- 
tions X i— > G t (k,y)XG t (k,y), corresponding to the identical operators (3.4), 
the reductions X i— > ^ f t [AT](y) are permutationally symmetric, and are not 
mixing only if the pure operations if they do not break this symmetry. 

The derived mixing property of the reduced stochastic dynamics t i— > 
p[X](t,v) obtained above for the corresponding statistical states 

p M [X}(t,v) = (r j M ,Q[X](t,v)r ] M )^Tr{Xg M (t,v)} (3.9) 
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gives an increase in the entropy 

a M (t,v) = -Tr{g M (t,v)\ng M (t,v)} 

for an ensemble of identical partides even under the condition of complete 
nondemolition observation. According to (3.9), the reduced density operators 
g M (t,v) for the system of M identical partides gives the probability density 

p M (t,v t ) = Tr{g M (t,v t )} = Ï> M {dv t ) /V™ {dv t ) 

of the output process Vt- Here = the probability measure on 

Tqo = Tqo x A°° defined by the Poisson measure (3.7). This means that the 
a posteriori density operator g M (t) is defined as a stochastic positive trace 
class operator normalized in the mean sense 

\\g M (t)h = jTr{g M (t,v)}pM(dv) = l. 
The density g M (t) satisfies the stochastic operator equation 



dg M (t) + i [H, g M (t)]dt = dn t 



i M 

^Y.G{k)g M {t)G\k)-g M {t) 



(3.10) 



which has a unique solution for every initial condition g M (0,v) = g^f given 
by the density operator g^f for the M-particle states ^[-X"] = Tr{X gff}. 
Let derive the differential equation (3.10) in the equivalent integral form 

p M [X](t) = p M [X(t)](t) + f í P M mX(t - r)](A)) - X(t - r)](r)dn t (dA) 

JO JA 

for an p™[X\ = (r} M ,Xr} M ), r] M G H M , where X(t) = U M (tYXU M (t), 

i M 

U M (t) = e~ ÍHMt / h , nX}(\) = — £ G\k, X)XG(k, X) . 



k=i 



Taking into account the fact that the stochastic integral (2.6) in this equation 
is simply a finite sum for every v G Tqo and t, one can write it as recursive 
operator equation 



r<t 



$ t [X](v)=X+ £ *rWX](r,\)-X](v), 

(r,y)ev 
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for a stochastic operation $t, defining the solutions to Eq. (3.10), as in (3.9), 
in terms of the composition (3.5). 

But such recurrence has the unique solution $ t (v) = ^t(yi) ° ^t{U2) ° • • 
defined for a v = (r, A) G as the chronological composition of the maps 
* t (r, à):Xh ^[X](r, À) if r < í and ^[X](r, A) = X if r > í. This solution 
can be found by the iterations 

$t(v)=I+ ]T $r(^)oA(r,A) 

(r,\)ev 

= I + E ( I+ E ° A( S , A)) o A(r, A) 

= . . . = £ A(zi) o . . . o A(z n ) , 

where A(y) = ^(y) — I, 3> s (i;) is the identical map I : X i— > X if s — ti and 
Eaçv t H z i) ° • • • ° A (^ t ) = %(yi) ° *t(j/2) o ... in terms of cr = . . . , ^ n }, 
2 = (r, A), si < . . . < s n , n < n t . 

Let us also write the nonlinear stochastic equation 

dtf(t) + i ^(í)]dí = eM(í) o - l)dn t (v) , 

where = #(A„ tH )/Tr{£(A„ tH )p}, 

i M i M 

£ ° = ^ E A)eGt(fc, A) , -E(A) = - £ G*(fc, X)G(k, A)) 
fc=i fc=i 

for the normalized density operator 

This describes the conditional expectations p^[X](t) = Tr{X^(t)} of the 
M-particle operators with respect to the output probability measure 

P M (dv) = P M (t,v)P^(dv), 

where p M {t,v) = Tr{g M {t, v)}. 
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4 Macroscopic and Continuous Reduction 
Limits 



We now consider the mean field approximation of the measurement apparatus 
fixing its total effect uk — 7 given by the mean number v of scattered bubbles 
(photons) per second and an interaction constant k coupling each bubble to 
a particle (atom) in the Hamiltonian (2.1). We look for the limits of the 
unitary and reduced evolutions (2.2) and (2.10) as v — > 00 and k — > such 
that 7 is a real constant. To perform these limits we need the expansions 



S(n) 
G(X) 



= I <g> 1 + i 



i^-R ® Q(n) - Q ) (£f (R ® Q(n)) 2 + . . . (4.1) 



of the scattering operator S(n) = exp{^ nR(E)Q(n)} and the reduced operator 
G(\) = fo(XI — k,R)/ fo(\) with respect to the coupling constant /t. The first 
term of the expansion for S(n) is disappearing in the r.h.s. of Eq.(2.2) while 
the second and third terms are appearing as the differentials of the operator- 
valued stochastic integrals 



h t [Q] 



Q(n r )dn r , n t [Q 2 



The corresponding terms 



JL 

-fo 



Ja -/o(A) 



dn r (dÀ) , n t 



fií 
Jo 



fo 



Q{n r ) 2 àn r . 



Ja /o(A) 



dn r (dÀ) 



in the right-hand side. of eq.(2.10) can also be written as the integrals 
n t [L] = Jq L(n r )dn r with vàlues in operator functions of r e 



nt[L\{r) 



n t (r) 

E L (n) 

n=l 



L(n) = l 0(n - 1} ® L ® 1 



(4.2) 



where L = [/(À)] is one of the multiplication operators 

L' = -[/í(A)//o(A)] , I/" = [fo(X)/ /o(A)] . 
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Here /ó(A) denotes the derivative dfo(X) = dfo(X)/dX and f"(X) denotes 
ddf(X). The stochastic integral n t [L](r) corresponding to the pointwise mul- 
tiplication L : f i— > lf of / e £ 2 (A) by a function / of the bubble coordinate 
A acts in £ as the multiplication operator 

nt(r) 

[n t [L](r)^](A) = £ Z(A n )^(A) = n t [Z](r, A)y?(A) . 

n=l 

Hence, the main terms on the right-hand side of Eq.(2.2) and (2.10) for 
k — > are given by the renormalised stochastic integrals 

X{t) = L{n r )dn r = —h t [L] (4.3) 
z/r Jo vt 

of the operator-valued stochastic functions L(t,r) = L{n t {r)) with respect 
to the numerical process n t (r) that has the Poisson probability distribution 
(2.5) on IV 

To pass to the large number limit v — > oo in (4.3) for an arbitrary 
operator L in L 2 (A), we need to use the quantum stochastic representa- 
tion [22] of the integral (4.3) in the Fock space T over L 2 (R + x A). The 
space JF can be defined as the L 2 (T)-space of all square integrable functions 
ip : T — > C, || ip\\ 2 = J r \(p(v)\ 2 X(dv) < oo of time ordered finite sequences 
v — (yi, . . . , y n ), y = (t, X) identified with subsets v C R+ x A of cardinality 

\v\ = 0,1,2, The measure X(dv) on the union T = S^=o^( n ) °f ^ ne 

disjoint subsets T(n) = {v G T : \v\ = n} is given as the sum 

oo 

A(A) = £A(T(n)nA) 

n=0 

of the product X(dv) = Y[ y&v dy of measures dy = dtdX on R + x A such that 
IM| 2 = E // \<p(y 1 ,...,y n )\ 2 l[dy i . 

n=Q J J0<ti<...<í„<oo i=1 

Let define N t [L] define the numerical integral in T by the action of the 
operator (4.2)in each Fock component T[r) = L 2 (A' T ') for all r e T: 

nt(r) 

\N t [L] V }{r) = L(n)<p(T) = n t [L]{r) V {r) . (4.4) 

n=l 
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Here ip(r) is the function </?(r, A) = (f(v) of A G A' T \ corresponding to a 
V? G L 2 (T) with a fixed time component of v — (r, A) G T. 

In order to obtain the initial probability measure P (di>) = 7r(dr)/ig°(dA) 
on Tqo induced by an initial Fock vector ip G L 2 (T), we need an isomorphic 
transformation of (4.4) 

N t [L] = N t [L] + MMtiL] + A\[Lf ]) + utf^Lf (4.5) 

which can be locally performed by a unitary transformation 

N t [L] = U}N t [L}U s , U s = exp{V^(4[/ ] - A a \fí])} 

for every t < s. Here and Al/ 1 ] are the creation and annihilation 

integrals of / G L 2 (A), p G L 2 (A)*, given by the operators 

'4[/Hl·l = E/( A Mf\v) 

A t [/t]J( v ) = / /(A)V(^Uy)dy 

in the Fock space L 2 (T), where v\y means the sequence v G T with deleted 
y = (r, X), r < t, and v U y means the sequence v G T with an additional 
element y ^ v. The characteristic functional of the stochastic operators n t [L] 
with respect to the initial state-vector /q° G £ and the Poisson probability 
measure (2.5) is now given simply by the vacuum expectation 

i (/ oO ) ei n t[ L] r )MdT) = {S e iN t[ L] S ) > 

jr°° 

where S^(v) — 1 if v — 0; otherwise, 5</,(f ) = 0. 

The corresponding representation /(í) = ^ A^ t [L] for (4.3) helps us imme- 
diately obtain the quantum large number limit 

hm - + N t [L] =fÍLf 1 

as the mean value l = (/ , Lf ) = fç,Lf of a single-bubble operator with 
respect to an initial wave packet / G £ 2 (A). This gives the following macro- 
scopic limit 

dip{t) + ^H^(t)àt = ^j(R ® q ï)ij;(t)dt 
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of the generalized Schròdinger equation (2.2) which turns out to be a nonsin- 
gular one with an additional potential — r jqoR corresponding to the mean mo- 
mentum q = (/ , Qfo) of a bubble in the initial state f . As one could expect, 
the mean field dynamics preserves the product structure ip(t,v) = r)(t)ip (v) 
of an initial product- vector if) = i] ® tp being trivial on the Fock component 
(po G T because of H — H <g> 1. But unexpectedly (compare with [19]) the 
macroscopic limit 

d X (t) + l - H X (t)dt = l - iq Rx(t)dt (4.6) 

of the reduction equation (2.10) corresponds to the same unitary dynamics 
r)(t) = U(t)r) = x(t) of the particle state- vector if x(0) = r], because of 

(fo,L'f ) = J / (A)VÓ(A)dA=i(to. 

The macroscopic limits for the M-particle system also give essentially the 
same continuous unitary evolutions in the large space H M ® T and in the 
reduced space TÍ M . To get this correspondence, one has only to replace the 
measure (2.5) on for Eq.(3.2) by the product measure tTq M on and 
for Eq.(3.10) by the induced measure (3.7) on T^, taking Mv instead of v in 
(4.5). This means that the mixing property of the reduced equation (3.10) 
vanishes in the mean field approximation for the bubble system. 

Let us now pay attention to the fluctuations with respect to the obtained 
large number limits. Such fluctuations might appear for k — jv — >0in the 
large time scale t ~ 1/k. We can get these fluctuations without rescaling the 
time t if we assume that qo = and k = so that we have to take into 

account also the /t 2 -terms in (4.1). 

It follows from the Fock space representation (4.5) that the quantum 
central limit 

hm -LN t [L]=A t [fU]+At[Lf ] 

exists for any single-bubble operator L with zero mean value (/o, Lf ) = 0. 
We first apply this central limit theorem to the right-hand side in (2.2): 

cty(í) + K ^(t)dt = ^(R® dut)^(í) . (4.7) 



vt 



Ió 



-fn 
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Here K = K <g> 1, K = i # + | ^ j ifo 2 /?, 

= Att/jQ] + 4 [Ç/o] = 2SMj[Q/ ] 

is the Fock space representation of the Wiener process u% with the dispersion 
a 2 = h 2 (dfo, dfo), defined by the quantum stochastic multiplication formula 
[21, 22] 

düdü = dA t \fÍQ]dAl[Qf ] = /ÍQQ/odí . 

The central limit equation (4.7) for the unitary evolution of the coupled 
system turns out to be a stochastic Schròdinger-Ito equation of diffusive 
type driven by the Wiener process Ut- The same conclusion obviously holds 
for the M-particle system driven by M independent Wiener processes üt(k) 
identical to ü t , k — 1, . . . , M: 

i / M \ 
dÍJ M (t) + K^màt = j.7[£R(k) ® àü t (k)y M (t) , 

where K M = ^H M + \ ^ E^li R{k)a 2 R{k), = K M ® ï. 

The application of the central limit theorem to the r.h.s. of the reduction 
equation (2.10) yields an essentially different type of the stochastic evolution, 

dx(í) + K x (t)dt = 1 R X (t)dv t , (4.8) 

originally derived in [21] by quantum calculus method. Here v t is a complex 
Wiener vector-process, with the Fock-space representation 

v t = A t [flL'] + A\[L'f ] = m\[(w + w )f ] + &4[(w - w )fo] (4.9) 

given by the complex osmotic velocity wo(X) = —d ln /o(A) of a single bubble, 
and the operator K is essentially the same as in (4.7). The linear stochastic 
equation (4.8) has a unique solution x{t-> v ) = T(t,v)r) for a given x(0, v) = 
i] G 7i which is not normalized ||x(í, f)|| 7^ for every Wiener trajectory 
t i— > Vt but is normalized in the mean square sense / \\x(t, ^)|| 2 Po(df ) = 1 
with respect to the Gaussian probability measure Po of v — {v t \t > 0}. The 
measure P is defined by the zero mean vàlues of v t and by the table 

dvdv = dA t [f^L')dA\[L'f ]=f^L'L'f dt 
dv*dv = dA t [f^]dA\[L'f ] = fÍL«L'f dt 
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of commuting multiplication 



dvdv* = dv*dv . 

But the reduction noise v t obtained in the Fock space representation does 
not commute with the real Wiener process ü t = ül in (4.7), 

dvdü = dA t [flL']dA\[QU}= flL'Qf ü dt 
düdv = dA t [flQ]dA\[L'f Q } = ÍÍQL'f dt 

if g(\) = d 2 ln/ (A) ^ 0, since 

[dv, dü] = fl[L', Q]f dt = * (/o, gf )dt . 

In the same way one can obtain the continuous reduction equation 

dg M (t) + (Kg M (t) + g M (t)KÏ)dt = 

( í ) E R{k)(r 2 Q M {t)R{k)dt + 1 {dw t Rg M {t) + g M (t)Rdw* t ) (4.10) 
k=i 

where R = J2if=i R{k) and the operator K = K M is the same as in the 
equation for the unitary evolution of the M-particle system coupled with the 
bubbles. The derived stochastic equation for the M-particle density operator 
g M (t,w) normalized in the mean is driven by the complex Wiener process 
Wt = vf 1 having the same multiplication table as the process \fMv t in (4.8). 
The diffusive type equation (4.10), (3.10) also has the mixing property. It 
was also derived in [5] by operational method. 



5 Conclusion 

The projection postulate for the reduction process follows as a result of con- 
ditioning by the results of a nondemolition measurement in the von Neumann 
Hamiltonian model [6] for the particle-meter interaction. This model gives 
also all unsharp reductions for the measurements with continuous spectrum. 

The singular interaction, corresponding to the collision model of scat- 
terings can be treated in quantum mechanics also in terms of Schròdinger 
equation in the generalized sense. The spontaneous localization [19, 20] of 
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the particle under the continual observation is explained as the result of the 
nondemolition countings of the scattered bubbles in a measurement appara- 
tus like cloud chamber. 

The spontaneous localization of a system of the identical partides (àtoms) 
in a bubble chamber (photodetector) is mixing due to the indistinguishability 
of the partides via the measurements of the bubbles. It is described in terms 
of the filtering equation for the density matrix of the partides, driven by 
the Poisson process of the total scatterings of the bubbles. This equation 
was derived in [5] within the operational approach to the quantum continual 
measurements. 

The macroscopic limit of the spontaneous localization equation is de- 
scribed in terms of the Schròdinger one, corresponding to the mean field 
approximation. The central limit of the spontaneous localization equation 
proves the diffusive reduction equation, which appeared in [17, 18] for some 
cases of R and H . Our theory shows the origin of the Wiener processes gener- 
ating the Ito diffusion equations for quantum states [23, 24], and also derives 
a diffusive mixing equation for the reduced density matrix of the system of 
partides under the continual observation. 

The quantum jump model, based on Poisson distributed impulsive mea- 
surements, is not new, see for example [25, 26], and it has been described on 
the statistical (a priori) level by the master equation which also results from 
the quantum diffusion models. However the treatment of these processes on 
the individual (a posteriori) level gains different types of quantum dynamics 
which can be described only in terms of the stochastic differential equations. 
The first such treatment of the counting measurements and the a posteriori 
dynamics was suggested in [13, 14] in terms of quantum stochastic calculus 
method [22]. The present work puts this method on the sòlid mathemati- 
cal basis by treatment of 5-function couplings with the classical Ito calculus 
applied to the quantum jumps and spontaneous localizations. 
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